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Sterile neutrinos with mass ~ 1 eV and order 10% mixing with active neutrinos have been proposed 
as a solution to anomalies in neutrino oscillation data, but are tightly constrained by cosmological 
limits. It was recently shown that these constraints are avoided if sterile neutrinos couple to a 
new MeV-scale gauge boson Al . However, even this scenario is restricted by structure formation 
constraints when yb-ntechahed collisional processes lead to efficient active-to-sterile neutrino conver¬ 
sion after neutrinos have decoupled. In view of this, we reevaluate in this paper the viability of 
sterile neutrinos with such “secret” interactions. We carefully dissect their evolution in the early 
Universe, including the various production channels and the expected modifications to large scale 
structure formation. We argue that there are two regions in parameter space — one at very small A 1 
coupling, one at relatively large Al coupling — where all constraints from big bang nucleosynthesis 
(BBN), cosmic microwave background (CMB), and large scale structure (LSS) data are satisfied. 
Interestingly, the large Al coupling region is precisely the region that was previously shown to have 
potentially important consequences for the small scale structure of dark matter halos if the Al boson 
couples also to the dark matter in the Universe. 

I. INTRODUCTION 

The possible existence of extra, “sterile”, neutrino 
species with masses at the eV scale and 0(10%) mix¬ 
ing with the Standard Model (SM) neutrinos is one of 
the most debated topics in neutrino physics today. Sev¬ 
eral anomalies in neutrino oscillation experiments m 
seem to point towards the existence of such particles, but 
null results from other experiments that did not observe 
a signal cast doubt on this hypothesis GHEE]- A multi¬ 
faceted experimental program is under way to clarify the 
issue and either detect, or conclusively rule out, eV-scale 
sterile neutrinos with large mixing angle. 

If the SM is indeed augmented with one or several such 
sterile neutrinos, but nothing else, some of the tightest 
constraints come from cosmological observations. In par¬ 
ticular, measurements of the effective number of relativis¬ 
tic particle species in the primordial plasma, N eS mm, 
disfavor the existence of an abundance of light or massless 
particles beyond the SM neutrinos and the photon in the 
early Universe. If sterile neutrinos are at the eV scale 
or above, they are also constrained by the distribution 
of large scale structure (LSS) in the Universe [14] which 
would be washed out due to efficient energy transport 
over large distances by free-streaming neutrinos. 

Cosmology, however, only constrains particle species 
that are abundantly produced in the early Universe. In 
two recent letters Of H], it was demonstrated that 
the production of sterile neutrinos can be suppressed 
until relatively late times if they are charged under a 
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new interaction. This idea has elicited interest in de¬ 
tailed model building and has interesting phenomenolog¬ 
ical consequences [MIS]- However, Mirizzi et al. [22] 
have recently pointed out that collisions mediated by the 
new interaction can result in significant late time produc¬ 
tion of sterile neutrinos and lead to tensions with CMB 
and LSS data on structure formation. There are, how¬ 
ever, several important caveats to this statement. In par¬ 
ticular, the bounds from Ref. [ 22 ] can be evaded if the 
sterile neutrinos either never recouple with active neutri¬ 
nos or remain collisional until matter-radiation equality. 
(We communicated on these caveats with the authors of 
Ref. [ 22 ]) w ho were aware of the first possibility but did 
not mention it as they found it less interesting in the 
context of previous literature. They mention the second 
possibility in the final version of their paper.) 


Our aim in the present paper is to understand in de¬ 
tail the role and impact of sterile neutrino collisions, and 
reevaluate if secretly interacting sterile neutrinos remain 
cosmologically viable. We begin in Sec. [Il] by review¬ 
ing the main features of self-interacting sterile neutrino 
scenarios. Then, in Sec. m we calculate the additional 
contribution to N eff at the BBN and CMB epochs. In 


Sec. IV we consider the impact on the large scale struc¬ 
ture in the Universe, focusing in particular on the sensi¬ 
tivity of Sloan Digital Sky Survey (SDSS) and Lyman-a 
data. We find that there are two regions of parameter 
space where sterile neutrinos with secret interactions are 
only weakly constrained. In Sec. [V] we discuss our con¬ 
clusions and summarize the results. 
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II. SECRET INTERACTIONS AND STERILE 
NEUTRINO PRODUCTION 

We assume that the Standard Model is augmented by 
a sterile neutrino v s with mass m J^| and with order 10% 
mixing with the SM neutrinos. We moreover assume the 
existence of a new secret t/(l) s gauge interaction, medi¬ 
ated by a vector boson A' of mass M at the MeV scale 
and coupling to sterile neutrinos through an interaction 
of the form 


Ant = ■ (1) 

Here, e s is th C/(l) s coupling constant and P L = (1 — 
y 5 )/2 is the left-handed chirality projection operator. We 
define the secret fine structure constant as a s = el/(Aw). 

This new interaction generates a large temperature- 
dependent potential [16] 
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( 2 ) 


for sterile neutrinos of energy E and sterile sector tem¬ 
perature T s . This potential leads to an in-medium mixing 
angle 0 m between active neutrinos v a and sterile neutri¬ 
nos v s , given by 


sin 2 2 9 m 


sin 2 20 O 

(cos 20 n + U eff ) + sin 29 0 


( 3 ) 


In the following, we will use a vacuum mixing angle 
0g ~ 0.1 and an active-sterile mass squared difference 
Am 2 ~ 1 eV 2 . As shown in mm, the secret interac¬ 
tions can suppress 0 m , and thus active to sterile neutrino 
oscillations, until after neutrino decoupling as long as 
\V eB \ » |Af7T 2 /(2i?)|. 

The new interaction also leads to collisions of sterile 
neutrinos. The collision rate for v s v s o v s v s scattering 
is given by 


Aoii n v a 


for T s » M 


where n v is the sterile neutrino density. The sterile neu¬ 
trino production rate T s and the final density depend on 
this collision rate. Two qualitatively different scenarios 
must be distinguished: 

Collisionless production: If the collision rate T co u is 
smaller than the Hubble rate H at all times, the ac¬ 
tive and sterile neutrinos can be taken to be oscillating 


Since v s is not a mass eigenstate, m s actually means the mass 
of the fourth, mostly sterile, mass eigenstate. 


without scattering pb If Am 2 /(2T Va ) H , v s are 

then produced only through oscillations, so that the final 
sterile neutrino number density is n„ ~ | sin 2 20 m n v , 

where n„ a = 3£(3) /{A'K 2 )g v is the density of one of 
active neutrino flavors and T v is the active neutrino tem- 
peratur. The final population of sterile neutrinos thus 
remains small, at most 0(10 2 ) of the active neutrino 
density, because of the small mixing angle. 

Collisional production: If T con exceeds the Hubble 
rate H , then sterile neutrinos cannot be treated as non- 
collisional [27]. In each collision, the sterile component of 
a v a -v s superposition changes its momentum, separates 
from the v a component, and continues to evolve indepen¬ 
dently. Subsequently, the active component again gen¬ 
erates a sterile component, which again gets scattered. 
This process continues for all neutrinos until eventually 
the phase space distributions of v a and v s have become 
identical. Thus, the fraction of v a converted to sterile 
neutrinos is not limited by the mixing angle, and all neu¬ 
trino flavors end up with equal number densities. 

The v a —> v s production rate in this case is T s ~ 
| sin 2 20 m • T coll [2Sj, where we can interpret the first 
factor as the average probability that an initially active 
neutrino is in its sterile state at the time of collision. 
The second factor gives the scattering rate that keeps it 
in the sterile state. We note that the production rate 
T s is proportional to n v and thus rapidly approaches its 
final value, 
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Note that, when T coU is much larger than the oscilla¬ 
tion frequency, using the average oscillation probability 
| sin 2 20 m is inappropriate, and in fact the production 
rate T s goes to zero in this case. Such a situation is, 
however, not realized for the parameter values explored 
in this work. 

In the following, we will look at both collisionless and 
collisional production of sterile neutrinos in more detail]^] 
with a special focus on the latter where more sterile neu¬ 
trinos may be produced. 


We ignore the SM matter potential and scattering experienced 
by active neutrinos because we will be interested in the regime 

where the secret interaction dominates over the SM interaction. 
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There is also the possibility that Mikheyev-Smirnov-Wolfenstein 
(MSW) type resonant effects, e.g., because of the sign-flip of 
the secret potential V e Q around T s ~ M, modify the v s pro¬ 
duction probability. In this work we treat all MSW transitions 
to be completely non-adiabatic and thus ignore them. A careful 
momentum-dependent treatment, which we defer to future work, 
is needed to accurately describe resonant conversion. 
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III. CONSTRAINTS ON A eff 


Cosmology is sensitive to the presence of relativistic de¬ 
grees of freedom through their contribution to the overall 
energy density. At early times the sterile sector was pre¬ 
sumably in equilibrium with the SM plasma, so that v s 
and Al were thermally populated. We assume that the 
sterile sector decouples from the SM sector well above 
the QCD scale and that oscillations remain suppressed 
until active neutrinos also decouple. Since the temper¬ 
ature T 7 of the SM sector drops more slowly than the 
sterile sector temperature T s when extra entropy is pro¬ 
duced during the QCD phase transition, T s at BBN is 
significantly smaller than T 7 . 

It is useful to track the ratio 

£ = Zsm ( 6 ) 

of the sterile sector temperature T s and the the temper¬ 
ature T v M of a standard neutrino. Before e + e annihi¬ 
lation, T® m = T 7 , while afterwards T® M = (4/ll) 1/,3 T 7 . 
Assuming comoving entropy is conserved, the ratio £ at 


BBN is 




£bbn = { 

( / 10.751 
1 V 106.7 J 

1 V 2-7/8 / 
i 

for M > 0.5 MeV 


| / 10.75 > 
I V 106.7 ^ 

i 3 

for M < 0.5 MeV 

J 

\ 0.649 

for M > 0.5 

MeV (case A) 

1 

[ 0.465 

for M < 0.5 

MeV (case B) 


(7) 

Here, the factor (10.75/106.7) 1,/3 gives the ratio of the 
sterile sector temperature to the active sector tempera¬ 
ture before A' decay, assuming that the two sectors have 
decoupled above the electroweak scale. It is based on 
counting the SM degrees of freedom that freeze out be¬ 
tween the electroweak and BBN epochs. 

A 1 is present in the Universe at the BBN epoch if 
M < 3T s | bbn ~ 0.5 MeV, and has decayed away if heav¬ 
ier. The factor (2 • 7/8 + 3) /(2 • 7/8) in the first row of 
Eq.|7f corresponds to the ratio of sterile sector degrees of 
freedonj^] before and after the decay of A 1 at T s ~ M/3. 

The extra radiation in the Universe is parameterized 
as N eS = (p v +p i/ A ')/p^, M , i.e., the energy density of all 
non-photon relativistic species, measured in units of the 
energy density of a SM neutrino species. The primordial 
population of v s and A' leads to N eS marginally larger 


Note that in complete models, for instance in scenarios including 
a dark Higgs sector to break the U(l) s symmetry, more degrees 
of freedom may need to be taken into account in the above equa¬ 
tions. 
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Figure 1. Possible cosmological histories of the active neutri¬ 
nos v a , the sterile neutrinos v s , and the sterile sector gauge 
bosons A 1 below the electroweak (EW) scale. Various possi¬ 
bilities, labeled as Al, A2 and Bl, B2, B3, are determined 
by the values of the Al mass M and the U( l) s fine structure 
constant a s and lead to testable predictions for A eff at both 
the BBN and CMB epochs. See text for details. 


than 3. For M > 0.5 MeV, 

Aeff,BBN(A) = N Va + (bBN(A) — 3.22 , (8A) 

at the BBN epoch. The first term, N v , on the right 
hand side accounts for the active neutrinos and is equal 
to 3.045. The second term includes the relativistic sterile 
sector particles, i.e., only v s if M > 0.5 MeV. If the A' 
bosons are lighter, i.e., M <C 0.5 MeV, they are present 
during BBN and contribute g A > = 3 degrees of freedom 
in the sterile sector, in addition to the g„ =2x7/8 
degrees of freedom of a sterile neutrino. Using the fact 
that also each active neutrino species has g„ =2x7/8 
degrees of freedom, we find 


^eff,BBN(B) — 


—Cbbn (B )-3.17, (8B) 


In Fig. Bl these two cases are summarized as BBN (A) 
and BBN(B), respectively. In either case, V efI BBN(A/B) 
remains consistent with the current BBN bound on extra 
radiation, AA r efrBnN = 0.661^45 (68%C.L.) [T2] , 

After BBN, the next important event is a possible se¬ 
cret recoupling of v a and v s . If the sterile neutrino pro¬ 
duction rate T s > H, a new hotter population of v s can 
be collisionally produced from v a , and they achieve ki¬ 
netic equilibrium with the primordially produced colder 
population of v s . Also, the Al can decay and heat up the 
sterile neutrinos. Both processes change the number and 
energy density of neutrinos, and N eS at CMB depends 
on the order in which they occur. 

In Fig. [2j we show the collisional v a -A v s production 
rate T s , normalized to the Hubble expansion rate 77, as 
a function of the photon temperature Z/: T s /H is sup¬ 
pressed at high temperatures (say, above GeV), where 
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sin 2 2 9 m is small due to the large V eS . Since F s oc T~ 3 in 
this regime (see Eqs. |2]), ^ and |5])) and H oc T 2 ,T S /H 
increases with T~ J as the temperature decreases. We de¬ 
fine the recoupling temperature T re as the temperature 
where T s /H > 1 for the first time since the primordial de¬ 
coupling of the active and sterile sectors above the QCD 
phase transition. When T s ~ M, the energy and temper¬ 
ature dependence of r s changes (see Eq. |5])), and when 
also V e ff drops below Am 2 /(2T S ) at T s < M, r s begins 
to drop again. The asymptotic behavior is E S /H oc T . 7 
at T s -C M and 9 m ~ 9 0 . There are then three possible 
sequences of events: 

1. No recoupling : For a sufficiently small interaction 
strength a s , the scattering rate T s always stays 
below the Hubble rate and there is no recoupling 
(solid black curve in Fig. [2]). 

If the interaction is stronger, a recoupling of v a and v s 
can happen either after or before A' decay: 

2. Recoupling after A' decay: If M > fewx 10 -2 MeV, 
the recoupling happens after A 1 have decayed (dot¬ 
ted blue curve in Fig. [2]). 

3. Recoupling before A' decay: If M < few x 
HE 2 MeV, the recoupling happens before A 1 have 
decayed (dashed red curve in Fig. [2]). 

In the second and third cases, there is also a secret de¬ 
coupling when r s /H again drops below one. If e 2 /M 2 < 
0(10 MeV -2 ), this decoupling happens while v s are still 
relativistic, i.e. T s > m s /3Q 

In the following, we discuss the three aforementioned 
cases in detail. 


1. No Recoupling 

In the no recoupling cases, labeled as Al and B1 in 
Fig. 0 the cosmological evolution after BBN is very 
straightforward. Vacuum oscillations convert a small 
fraction ~ \ sin 2 2 9 0 ~ 0.01 of active neutrinos into ster¬ 
ile neutrinos (and vice versa), but this has negligible 
impact on the cosmological observables. Therefore, the 
temperature ratio £ at CMB can be derived from the sep¬ 
arate conservation of entropy in the active neutrino sec¬ 
tor and in the sterile neutrino sector. It is independent of 


5 In this paper, we will always assume this to be the case since we 
will find that the parameter region with e 2 /M 2 >0(10 MeV -2 ) 
is already disfavored by the requirement that active neutri¬ 
nos should free stream sufficiently early m (see Secs. |IV| and 
0 - If u s and u a are still coupled when the v 3 become non- 
relativistic, the mostly sterile mass eigenstate v 4 will undergo a 
non-relativistic freeze-out and partly annihilate to pairs of mostly 
active neutrinos. Similarly, there is the possibility that the A de¬ 
cay after the decoupling, but this does not happen for the range 
of parameters we will discuss here. 



Figure 2. Evolution of the collisional v a —» v s production rate 
r s , normalized by the Hubble rate H, versus the photon tem¬ 
perature T~, for different representative choices of the secret 
gauge boson mass M and the secret fine structure constant 
a 3 . When T s /H > 1 , collisional production of v s from the 
thermal bath of v a is effective. The solid black curve shows a 
case where this never happens. The shoulder around T 7 ~ M 
is where the A 1 decay away. The dotted blue and dashed red 
curves correspond to recoupling after and before A' decay, 
respectively. 


when the A' decay (provided that it happens before the 
CMB epoch and approximately in chemical equilibrium.) 
That is, (cmba/b — £bbn(aj = 0.649. 

lV eff cMB can be estimated in analogy to Eq. (8A|. 
For the assumed sterile neutrino mass ~ 1 eV, the v s 
contribution to the relativistic energy density has to be 
weighted by an extra factor because they are already 
semi-relativistic at the CMB epoch, where the photon 
temperature is T 7 ~ 0.30 eV and the kinetic temperature 

of the sterile sector is T s = £cmb ■ (d/ll) 1 ^ 3 ^ ~ 0.14 eV. 
As in [22] , we assume that the extra weight factor is char¬ 
acterized by the pressure P. (See Appendix [A] for the 
definition and calculations of the kinetic temperature and 
the pressure P used here.) We thus obtain 


^eff,CMB — N Va -)-— 5 -— • £cMB — 3.13 . (9) 

ma=0 CMB 

It is worth noting that the CMB temperature spectrum 
does not exactly measure the value of V effjCMn . Instead, 
the observed spectrum depends on the evolution of the 
energy density in relativistic degrees of freedom between 
the epoch of matter-radiation equality (T 7 eq ~ 0.7 eV) 
and recombination (T 7iCMB ~ 0.30 eV) HI]. Therefore, 
the value of AT eff CMB measured from the CMB temper¬ 
ature power spectrum lies between the values of N eS at 
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T 7 c mb an< i ^y.eq- The latter value, which we will denote 
by -/V effjeq , is given by 


N „ = N + 

1 v eff,eq v a ' 


m s = 1eV 


p. 


m= 0 


£eq — 3.18 . 


( 10 ) 


the time of matter-radiation equality: 

N eff , eq = 3.09. (14) 

Both values are in agreement with the Planck bound m- 


Both of these values agree with the bound from the 2015 

Planck data release, N eS = 3.15 ± 0.23 (68% C.L.) PJ. 3. Recoupling before A' decay 


2. Recoupling after A' decay 


The cases of recoupling after A' decay are labeled as A2 
and B2 in Fig. [I] In both cases, entropy conservation in 
the sterile sector before recoupling leads to a temperature 
ratio just after A' decay of ~ £bbN(a) = 0.649, which 
in turn implies 

= 3.045 + & = 3.22. (11) 

Here, we have assumed that during A' decay chemical 
equilibrium holds in the sterile sector. 

After recoupling, efficient neutrino oscillations and 
collisions lead to equilibration of the number densities 
and energy densities of all active and sterile neutrino 
species. Nevertheless, since number-changing interac¬ 
tions are strongly suppressed at T s -C M, this recou¬ 
pling cannot change the total (active + sterile) neutrino 
number density and energy density beyond what is ne¬ 
cessitated by cosmological expansion. The kinetic tem¬ 
perature T v re shared by all neutrinos after recoupling is 
then given by 


„ 3.045 • (Q 4 + 1 • Tg re 0 

~ 3.045 • C® 3 + 1 • T 3 rei0 


0.97 T™, 


( 12 ) 


where T^ M is the active neutrino temperature just prior 
to recoupling (which is at its SM value) and T s ie 0 de¬ 
notes the sterile sector temperature just prior to recou¬ 
pling. 

Eventually, the mostly sterile eV-scale mass eigenstate 
decouples from the light mass eigenstates and becomes 
semi-relativistic at the CMB time. Its kinetic temper¬ 
ature at this epoch is T Sj cmb — 0.13 eV. The effective 
number of relativistic species at the CMB epoch is given 

by 


^eff,CMB — -W- 


eff,M 


3 1 Pm s =l eV 

4 + 4 P m _ 0 


CMB 


~ 2.51 + 0.39 ~ 2.90. 


(13) 


Note that this is smaller than the SM value 3.045. This 
happens because part of the energy of active neutrinos 
has been transferred to the mostly sterile mass eigenstate 
i/ 4 , whose kinetic energy gets redshifted away more effi¬ 
ciently after it becomes non-relativistic. Ref. [22] also 
found N e g < 3 for this scenario. Similarly, we obtain for 


The last possibility, labeled as case B3 in Fig. [l| is that 
recoupling happens before A' decay. In this case, all neu¬ 
trinos, together with A 1 , reach a common chemical equi¬ 
librium, which lasts until most of the A' particles have 
decayed. During the formation of chemical equilibrium, 
the total energy is conserved while entropy increases. En¬ 
ergy conservation allows us to calculate the temperature 
Tj, re of the active + sterile neutrino sector immediately 
after recoupling: 


(3 Qi/ a + 9v s 4“ 9 a' ) Tjqre 

= [3 9v a + (flV, + 9a') £bBN(b)] (T(?re ) 4 , (15) 

where is again the active neutrino temperature just 
prior to recoupling. Plugging in numbers for the effective 
numbers of degrees of freedom g u , g v , g A < and using 
(bbn(b) = 0.465, we obtain 

7^3 = 0.8617^. (16) 

Later, the A 1 decay and the thermal bath of neutrinos 
gets reheated by a factor [(3 g v + g Vs + 5+)/(3<?!y a + 
g u j] 1 / 3 ~ 1.125. The effective number of relativistic 
species after A' decay is then 

N effM ~ + gt/ ° (1.125 • 0.861) 4 ~ 3.568. (17) 

9o a 

The next steps are the decoupling of sterile neutrinos and 
active neutrinos, and then the freeze-out of sterile neu¬ 
trino self-interactions. Since the number densities and 
energy densities of the different species do not change 
during these decouplings, the total effective number of 
relativistic degrees of freedom at the CMB epoch is given 
by 


^cff,CMB — -W- 


3 1 Pm s =le V 

^'4 + 4^r 


CMB 


3.21, 

(18) 


where again the pressure characterizes the contribution 
of the semi-relativistic iq to the radiation density in the 
Universe. Similarly, at matter-radiation equality we have 

iV effjeq ~3.43. (19) 

This number is still within the 2er error of the Planck 
bound [TS] . 
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IV. STRUCTURE FORMATION 


Besides the constraints on extra radiation measured 
by N e fi, CMB data also prefers that most of the active 
(massless) neutrinos start to free-stream before redshift 
z ~ 10 J [29|13Qj. On the other hand, matter power spec¬ 
trum observations forbid these free-streaming degrees of 
freedom from carrying so much energy as to suppress 
small scale structures m- Therefore measurements of 
the matter power spectrum put the most stringent up¬ 
per bound on the mass of all fully thermalized neutrino 
species: ')h m v ~ 0.2-0.7 eV (95% C.L.) ;T3]. This con¬ 
cern [22] excludes a large proportion of the parameter 
region for self-interacting sterile neutrinos considered in 
m■ However, like the constraint on N eS discussed in 
Sec. |III[ it is avoided if never exceeds H after the 
epoch when V e g drops below the oscillation frequency 
(cases Al and B1 in Fig. [Tj). 

Interestingly, structure formation constraints are also 
significantly relaxed when the [/(l) s gauge coupling e s is 
large and/or the gauge boson mass M is small. In this 
case, sterile neutrinos, although produced abundantly 
through collisional production (see Sec. [II]), cannot free- 
stream until late times, long after matter-radiation equal¬ 
ity. Thus, their influence on structure formation is sig¬ 
nificantly reduced. We will now discuss this observation 
in more detail. 

After the active and sterile neutrinos have equilibrated 
through W-mediated collisions, they should be treated 
as an incoherent mixture of the four mass eigenstates v i 
[i = 1.. .4). The reason is that for m 4 ~ 1 eV, their 
oscillation time scales are much smaller than both the 
Hubble time and the time interval between scatterings. 
For simplicity, assume that only the mostly sterile mass 
eigenstate zq is massive with mass m 4 ~ 1 eV, and that 
it only mixes appreciably with one of the mostly active 
mass eigenstates, say zq: 

zq ~ sin 9 0 zq + cos 6 0 zq . (20) 

We take the vacuum mixing angle to be 9 0 ~ 0.1 and 
we take into account that matter effects are negligible 
at temperatures relevant for structure formation (after 
matter-radiation equality). Since it is the flavor eigen¬ 
state zq that is charged under U(l) s , the mass eigenstates 
zq and zq interact with relative rates sin 2 9 0 and cos 2 9 0 , 
respectively, while zq and zq essentially free-stream. 

To study the influence of the secret interaction on 
structure formation, we estimate the mean comoving 
distance A s that each zq can travel in the early Uni¬ 
verse. Since neutrinos can transport energy efficiently 
over scales smaller than A s , the matter power spectrum 
will be suppressed on these scales. As long as neutrinos 
are collisional, they do not free stream, but diffuse over 


scales of order [32] 


(A™ 11 ) 2 ^ 


dt 


(Vs ) 2 

a 2 (t) 




( 21 ) 


where a(t) is the scale factor of the Universe, t dec is the 
time at which sterile neutrino self-interactions decouple, 


4 2 a 

( crv) a ~ (v s ) £s cos —AU (T + m s ) 2 

X ' x \m 2 + t 2 ) 2K ’ 


( 22 ) 


is the thermally averaged interaction cross section of 
the mostly sterile mass eigenstate zq, estimated here by 
naive dimensional analysis, and n s is the number density 
of sterile neutrinos. For simplicity, we take the kinetic 
temperature T s of the sterile sector equal to the active 
neutrino temperature in this section, i.e. T s = T^ M = 
(4/ll) 1//3 T 7 oc a _1 (f), as long as zq are relativistic. If 
sterile neutrinos become non-relativistic (T s < m s ) while 
they are still strongly self-coupled, the kinetic temper¬ 
ature of the sterile sector scales as T s oc a~ 2 (t) until 
T s drops below T sdec . After that, the sterile neutrino 
momenta are simply redshifted proportional to a 1 (t). 
This implies in particular that, at T s m s , we have 
n s ~ T 3 , while after jq become non-relativistic, but are 
still strongly coupled, this changes to n s oc T 3 ^ 2 . The 
computation of the average velocity (v s ) of iq entering 
eq. (221 is discussed in Appendix [Aj 


The decoupling temperature T s dec and the correspond¬ 
ing time t^ ec are defined by the condition that the sterile 
neutrino interaction rate is just equal to the Hubble rate: 


(<™> a 


= H(t dec ). 


(23) 


.dec 


sterile neutrinos start to free stream. The total 

dec 


After t 

comoving distance that a iq travels between the time f' 
and the present epoch t° is [32] 


A? = 


f dt Ml 

(dec a[t) 


The overall damping scale is then given by 

A 2 = (Af 1 ) 2 + (A?) 2 


(24) 


(25) 


At scales larger than A s , structure formation is unaffected 
by the existence of sterile neutrinos, while at smaller 
scales, structures are washed out. 

As a numerical example, for M = 0.1 MeV, e s = 0.1, 
we find A™ 11 ~ 29 Mpc /h, A^ s ~ 68 Mpc /h and thus 


\ s ~ 74 Mpc /h, 
corresponding to a wave number of 

k s = 2tt/X s ~ 0.085 h/ Mpc. 


(26) 


(27) 
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Figure 3. (a) Three-dimensional matter power spectrum Pm (|fc|) derived from the SDSS Luminous Red Galaxy (LRG) Sam¬ 
ple [ 33 ] and (b) one-dimensional flux power spectrum A 2 (k) = kP F (k)/n of Lyman-a photons at various redshifts, compared 
to the qualitative predictions of sterile neutrino models with (red dashed curves) and without (green dotted curves) self¬ 


interactions. Note that the relation between Pm (|fc|) and P F (k) is non-linear, see e.g. [34]. The assumed self-interaction 
parameters are e s = 0.1, M = 0.1 MeV, and the assumed sterile neutrino mass is 1 eV. The data points and the SM prediction 
(solid green curves) are taken from [33] and and from [35], respectively. The predictions including sterile neutrinos are obtained 
by multiplying the SM predictions by the fc-dependent suppression profile from Fig. 7 of EH, shifted such that the onset of the 
suppression is at our calculated damping scale k s (Eq. ( |27| ) ; see text for details), and scaled such that the maximum suppression 
is given by Eq. (291 for panel (a) and by (301 for panel (b). Note that the error bars shown here are statistical only, and large 
systematic uncertainties, especially at small scales (large k ) should be kept in mind. 


This should be compared to the free streaming scale of a 
decoupled sterile neutrino with a mass < 1 eV, 


no self-int. 

'S 


0.018 



ft./Mpc. 


(28) 


This factor of ~ 5 decrease in the free streaming scale 
compared to a conventional sterile neutrinos without self¬ 
interactions implies that data on large scale structure 
(LSS) and baryon acoustic oscillations (BAO) will be in 
much better agreement with our model than with ster¬ 
ile neutrino models that do not feature self-interactions. 
The strongest constraints will come from data probing 
very small scales, in particular Lyman-a forests. 


Even at scales k > k s , the suppression of the matter 
power spectrum Pm (| k |) does not set in abruptly, but in¬ 
creases gradually. For non-interacting sterile neutrinos, 
numerical simulations show that the suppression satu¬ 
rates at k ~ 50 k s . At even smaller scales (even larger 
k), the deviation from the prediction of standard cosmol¬ 
ogy is hqimi 


SP M (\k\) 

P M (\k \) 


( 29 ) 


in the linear structure formation regime. Here, /„ = 
3m s C(3)/(27r 2 )T s 3 (t 0 )x87rG/(3iT 2 (t 0 ))/H m ~ 0.07 is the 
ratio of the sterile neutrino mass density to the total 
mass density fl m ~ 0.3 today. In the regime of non¬ 
linear structure formation, 5P M (\k\)/P M (\k\) is some¬ 
what larger EH Eg, but N-body simulations show that 
it decreases again at scales k > few h/Mpc [34] [37]. 

It is, however, difficult to directly measure P M (\k\) 
at these nonlinear scales. The most sensitive data 
sets are Lyman-a forests, from which the 1-dimensional 
flux power spectrum Pp(k) of Lyman-a photons can 
be extracted. Translating Pp{k) into a measurement of 
P M {\k\) requires a determination of the bias b(k), which 
is obtained from numerical simulations of structure for¬ 
mation that include the dynamics of the gas clouds in 
which Lyman-a photons from distant quasars are ab¬ 
sorbed. For SM neutrinos, such simulations have been 
performed for instance in [341 , and we can estimate from 
Fig. 13 of that paper that the maximal suppression of 
P F {k) is of order 


SP F (k) 
P F (k ) 


- 0.1 x 


y 1 eV 


(30) 


where m „ is the sum of all neutrino masses. This esti¬ 
mate is crude but conservative, and ignores the fact that 





















the maximal suppression is actually smaller at lower red- 
shifts. The suppression of P F (k) described by Eq. (30) 
is smaller than the suppression of P M (|fc|) from Eq. (29) 
because of the nonlinear fc-dependent relation between 
the matter power spectrum and the flux power spectrum 
(see for instance [55] . especially Fig. 16 in that paper). 
Since no dedicated simulations are available for our self¬ 
interacting sterile neutrino model, we will assume in the 
following that 5P F (k)/P F (k ) saturates at the value given 
by Eq. (301 even when ^ m u is dominated by the ster¬ 
ile neutrino mass mn s . This amounts to assuming that 
the impact of secretly interacting sterile neutrinos on 
these small scales is qualitatively similar to that of ac¬ 
tive neutrinos. A more detailed treatment requires a 
dedicated simulation including these secretly interacting 
sterile neutrinos. Note that neutrino free-streaming after 
CMB decoupling may lead to less suppression than de¬ 
scribed in Eqs. (291 and (301 because perturbation modes 
well within the horizon have already grown significantly 
by that time. We will not include this effect in the fol¬ 
lowing discussion to remain conservative. 


From Fig. [3] (a), we observe that the suppression of 
the matter power spectrum at scales < 0.2 h /Mpc due to 
self-interacting sterile neutrinos is completely negligible, 
while a fully thermalized non-interacting sterile neutrino 
with the same mass leads to a clear suppression already 
at these scales. This implies that self-interacting ster¬ 
ile neutrinos with the parameters chosen here are not 
constrained by data on linear structure formation. Go¬ 
ing to smaller scales or larger k (Fig. [3] (b)), where non¬ 
linear effects become relevant, we see that both the sterile 
neutrino model with self-interactions and the one with¬ 
out lead to suppression, but the amount of suppression 
is reduced in the self-interacting case. It was shown in 
Ref. [35] that the data disfavors suppression larger than 
10% at k = 10/i/Mpc. Self-interacting sterile neutri¬ 
nos at the benchmark point shown in Fig. [3] appear to 
be marginally consistent with this constraint. It should 
be kept in mind, however, that our predictions are only 
qualitative. Therefore, only a detailed fit using simula¬ 
tions of non-linear structure formation that include ster¬ 
ile neutrino self-interactions could provide a conclusive 
assessment of the viability of such a scenario. 


We show the qualitative impact of self-interacting 
sterile neutrinos on large scale structure in Fig. [3] 
Panel (a) compares theoretical predictions in models 
with and without sterile neutrinos to data on the 
three-dimensional matter power spectrum P M (|fc|) from 
the Sloan Digital Sky Survey (SDSS) Luminous Red 
Galaxy (LRG) catalog [33] . Panel (b) compares to one¬ 


dimensional flux power spectra A 2 {k) = k P F (k)/TT from 
Lyman-a forest data [35]. SDSS-LRG data corresponds 
to a mean redshift of z ~ 0.35, while Lyman-a data is 
split up according to redshift and reaches up to z ~ 5.4. 
Note that the data in (351 is presented as a function of the 
wave number k v in velocity space, measured in units of 
sec/km. The conversion to the wave number k in coordi¬ 
nate space, measured in units of h/ Mpc, is done accord¬ 
ing to the formula k = k v H(z)/(l + z), where H{z) is the 
Hubble rate at redshift z. The theoretical predictions for 
the SM with vanishing neutrino mass (solid green curves 
in Fig. [3]) are taken from [33] and [35], respectively. Our 
(qualitative) predictions for sterile neutrino models with 
and without self-interactions are obtained in the follow¬ 
ing way: we start from the numerical prediction for the 
neutrino-induced suppression of the matter power spec¬ 
trum from Ref. [3Tj . In particular, we use the curve cor¬ 
responding to /„ = 0.07 from Fig. 7 in that paper. We 
then shift this curve such that the onset of the suppres¬ 
sion coincides with our calculated damping scale k s , and 
we rescale it such that the maximal suppression is —8 f v 
in Fig. [3] (a) (linear regime) and 10% in Fig. [ 3 ] (b) (non¬ 
linear regime), see Eqs. ([29]) and (30). We then multiply 


with the SM prediction to obtain the dotted green curves 
for sterile neutrinos without self-interactions and the red 
dashed curves for sterile neutrinos with self-interactions 
in Fig. [3j We use e s = 0.1, M = 0.1 MeV and m s = 1 eV. 
Since we neglect a possible upturn of the power spectrum 
at k > 1 h/Mpc [37], our estimates are very conservative. 


Let us finally discuss how the cosmological effects 
of the three active neutrinos are modified in the self¬ 
interacting sterile neutrino scenario. The dynamics of 
the mass eigenstates v 2 and ^ 3 , which we assume not to 
mix with zz 4 , is the same as in standard cosmology: they 
start to free stream at redshift z 10 °. v 1} however, 
starts to free stream later than a non-interacting neu¬ 
trino, but earlier than z/ 4 . The free-streaming condition 
for is, in analogy to Eq. ([23]) , 


(T S M ) 3 .( T S M )2( J \ s ^ Qq < H _ (31) 


As shown in Ref. [59], free-streaming of active neutrinos 
before redshift z ~ 10° is required to sufficiently suppress 
the acoustic peaks in the CMB power spectrum. The 
change from three to two truly free streaming neutrino 
species in our model will lead to minor modifications of 
the CMB power spectrum, but the analysis from :29j sug¬ 
gests that these are unlikely to spoil the fit to CMB data, 
in particular since they may be compensated by changes 
in the best fit values of other cosmological parameters. 


Note that also Planck CMB data alone, without includ¬ 
ing data on large scale structure observations, imposes an 
upper limit on the mass of sterile neutrinos, which, for a 
fully thermalized species is m s < 0.5 eV at 95% C.L. T3] , 
A much weaker bound is expected if self- interactions 
among sterile neutrinos are so strong that they remain 
collisional until after the CMB epoch. In this case, the 
early Integrated Sachs-Wolfe (ISW) effect induced by v s 
perturbations at low multipole order (50 < l < 200) will 
be suppressed m- Thus, the main effect on the CMB 
will come from the shift of matter-radiation equality, to 
which the sensitivity is, however, much weaker. 
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V. DISCUSSION AND CONCLUSIONS 


As we have seen in the previous section, there are two 
main scenarios in which self-interacting sterile neutrinos 
do not run into conflict with cosmological data: 

(i) The zq production rate T s drops below the Hub¬ 
ble expansion rate H before the effective potential |14fj| 
drops below the oscillation frequency |Am /(2U)| and 
the dynamic suppression of active-sterile mixing due to 
V eff ends. In this case, sterile neutrinos are not produced 
in significant numbers in the early Universe and hence 
cosmology is not sensitive to their existence. An expla¬ 
nation of small scale structure anomalies as advocated 
in [16] is, however, not possible in this scenario]^] In par¬ 
ticular, even if the new interaction also couples to dark 
matter as proposed in [TB] . it is too weak to have phe¬ 
nomenological consequences. 

This disadvantage can be avoided if more than one self¬ 
interacting sterile neutrino exists. Consider for example, 
a model with three mostly sterile neutrino mass eigen¬ 
states zq, zq, zq. Let zq and zq have a relatively large 
mixing 9 0 ~ 0.1 with active neutrinos, as motivated for 
instance by the short baseline oscillation anomalies. Let 
their coupling to the A! gauge boson be e^ 4 ’ 5 ^ — 10 -5 , 
large enough to dynamically suppress their mixing with 
the mostly active mass eigenstates until after BBN, but 
small enough to prevent their equilibration afterwards. 
On the other hand, let zq have a vanishing mixing with 
zq 2 , 3 ) but a larger secret gauge coupling ~ 0.1. Due 
to its small mixing, it is never produced through oscil¬ 
lations. However, its primordial population—the relic 
density produced before the visible and sterile sectors 
decoupled in the very early Universe—still acts as a ther¬ 
mal bath to which the dark matter may be strongly cou¬ 
pled, thus potentially solving the missing satellites prob¬ 
lem 

(ii) The self-interactions are so strong that sterile neu¬ 
trinos remain collisional at least until matter-radiation 
equality. In this scenario, sterile neutrinos are prod uced 
when |U eff | < |Ato 2 /(2£’)|. However, as shown in Sec. 


Ill 


the effective number of relativistic degrees of freedom in 
the Universe, N e g, remains close to 3 because equilibra¬ 
tion between active and sterile neutrinos happens after 
neutrinos have decoupled from the photon bath. More¬ 
over, as argued in Sec. |IV| the impact of self-interacting 
sterile neutrinos on structure formation is much smaller 
in this scenario than the impact of conventional non¬ 
interacting sterile neutrino because they cannot trans¬ 
port energy efficiently over large distances due to their 
reduced free-streaming. Structure formation constraints 
could be further relaxed in models containing, besides 
an eV-scale mass eigenstate zq, one or more additional 


Note that recent simulations of cosmological structure formation 
suggest that these anomalies may be resolved once baryons are 
included in the simulations |391140| . 


mostly sterile states with much lower masses [25]. It is in¬ 
triguing that the parameter region corresponding to sce¬ 
nario (ii) contains the region where small scale structure 
anomalies can be explained, as shown in |16| . 

We summarize these results in Fig. [4] The yellow cross- 
hatched region on the right is ruled out because active 
and sterile neutrinos come into thermal equilibrium be¬ 
fore the active neutrino decoupling from the SM plasma. 
In the lower part of this region, this happens simply be¬ 
cause is negligibly small. In the upper part, U e ff is 
large, but also r s is large so that collisional production of 
sterile neutrinos is efficient in spite of the suppressed in¬ 
medium mixing 9 m . This leads to constraints from N e g 
and from the light element abundances in BBN [2T|. In 
the blue vertically hatched region, sterile neutrinos are 
produced after zq decoupling, so that CMB constraints 
on N e g remain satisfied. However, sterile neutrinos free- 
stream early on in this region and violate the CMB and 
structure formation constraints on their mass. This mass 
constraint can be considerably relaxed if the sterile neu¬ 
trinos remain collisional until after the CMB epoch at 
T 1 ~ 0.3 eV. This defines the upper edge of the blue 
hatched region. In the red shaded region in the upper 
left corner, the secret interaction is too strong and zq 
free streams too late. CMB data requires that active 
neutrinos free stream early enough, and thus strongly 
disfavors this region. Two white regions remain allowed: 
Scenario (z) with weak self-interactions, corresponds to 
the wedge-shaped white region in the lower part of the 
plot. Scenario (ii), with strong self-interactions, is real¬ 
ized in the thin white band between the blue vertically 
hatched region and the red shaded region. As explained 
above, whether or not this white band is allowed depends 
strongly on systematic uncertainties at Lyman-a scales 
and on the possible existence of additional states with 
masses <C 1 eV. 

There are several important caveats and limitations 
to the above analysis. First, we have only worked with 
thermal averages for the parameters characterizing each 
particle species, such as energy, velocity, pressure, etc. 
To obtain more accurate predictions, it would be nec¬ 
essary to solve momentum-dependent quantum kinetic 
equations. This would be in particular interesting in the 
temperature regions where V e g changes sign and where 
V eS ~ Am 2 / (2T V ) x cos20 o . We expect that our mod¬ 
eling of flavor conversions in this region as fully non- 
adiabatic transitions is accurate, but this assumption re¬ 
mains to be checked explicitly. Moreover, the impact of 
self-interacting sterile neutrinos on non-linear structure 
formation at the smallest scales probed by Lyman-a data 
should be calculated more carefully. Improving these is¬ 
sues is left for future work. 

In conclusion, we have argued in this paper that self¬ 
interacting sterile neutrinos remain a cosmologically vi¬ 
able extension of the Standard Model. As long as the self¬ 
interaction dynamically suppresses sterile neutrino pro¬ 
duction until neutrinos decouple from the photon bath, 
the abundance produced afterwards is not in conflict with 
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Figure 4. Schematic illustration of the parameter space for 
eV-scale sterile neutrinos coupled to a new “secret” gauge 
boson with mass M and a secret fine structure constant a B . 
The vacuum mixing angle between active and sterile neutri¬ 
nos was taken to be 9 0 = 0.1. The white region in the lower 
half of the plot is allowed by all constraints, while the narrow 
white band in the upper left part satisfies all constraints ex¬ 
cept possibly large scale structure (LSS) limits from Lyman-a 
data at the smallest scales. The red stars show representative 
models in scenarios (i) and (ii). The colored regions are ex¬ 
cluded, either by LSS observations (blue vertically hatched), 
by the requirement that active neutrinos should free stream 
early enough (red shaded), or by a combination of CMB and 
BBN data (yellow cross-hatched). 


constraints on N e g. Moreover, if the self-interaction is ei¬ 
ther weak enough for scattering to be negligible after the 
dynamic mixing suppression is lifted, or strong enough to 
delay free streaming of sterile neutrinos until sufficiently 
late times, also structure formation constraints can be 
avoided or significantly relaxed. 
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Appendix A: Kinetic temperature and pressure 


In the following, we give more details on the momen¬ 
tum distribution function /(p, t) of sterile neutrinos v s 
after they have decoupled from all other particle species. 


/(p, t) is essential in the calculation of the pressure P in 


Sec. Ill and the average velocity (v s ) in Sec. IV 


Even when v s are decoupled from other particles, they 
may still couple to themselves via strong self-interactions. 
If the self-interaction conserves the number of particles, 
such as last's -o- v s v s , it only maintains kinetic equilib¬ 
rium, but not chemical equilibrium. Number conserva¬ 
tion and entropy maximization force the v s momentum 
distribution function in kinetic equilibrium to take the 
form 


e [E{p)-n s (t)\/T s (t) + 1 ’ ( A1 ) 

where T s (t ) is defined as the kinetic temperature , [i s (t) is 
the chemical potential, and E(p) = (p 2 + m 2 ) 1 ^ 2 . Here 
and in the following, we use the definition p = |p|. Since 
we are interested in the evolution at relatively late times, 
when the sterile neutrino density is low compared to the 
density of a degenerate fermion gas and thus (f(p,t)) <C 
1 , the classical approximation 

/(p,i)~e _[B(p)_ ' 1 * (t)]/T * (t) (A2) 


is adequate. 

Our goal is to solve for the functions T s (t) and p s (t) 
with the initial condition of a relativistic thermal en¬ 
semble of sterile neutrinos. This means that, initially, 
T s = Tj rn s and p s = 0 at a = a t . Note that the 
sterile neutrino mass will lead to a non-zero p s soon af¬ 
ter neutrinos go out of chemical equilibrium. Although it 
is difficult to analytically solve the corresponding Boltz¬ 
mann equation, there are two conditions that can be used 
to numerically obtain T s and p s as functions of the scale 
factor a(t). One is number conservation. The other is 
entropy conservation, which holds approximately for ki¬ 
netic equilibrium in the classical limit [iH] . 

The number density is 


>3 

a p 


■ m = J W? 


and the classical entropy density is defined as 


73 

a p 


(A3) 


3 (*)=/tVA3 /Cp> *) [! — ln *)] - ( A4 ) 

J y Zttj 


It is straightforward to obtain the asymptotic solu¬ 
tions [TE> 


and 


T s (t ) oc 



for T s m s 
for T s <C m s 


p s (t) oc 


a{t) 

const 


for T s m s 
for T s -C m s 


(A5) 


(A6) 
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In the transition region T s ~ m s , the solution needs to 
be obtained numerically. The result is plotted in Fig. [5} 
Finally, we comment on the calculation of the pressure 
P and the average velocity (v a ) of sterile neutrinos. The 
pressure is given by [T6j 


r 73 : 

p = I d p p 


= -T e pJT ’ [ d3p - — e - E{ - p)/T ‘ 

J (2 tt) 3 * 3 dp 


(A7) 


and the average velocity is 


{v s ) ^ 


1 

N 


d 3 p p 

(27r) 3 E(jp) 




(AS) 


Figure 5. Kinetic temperature T g and chemical potential 
p, s of sterile neutrinos, as functions of the scale factor a(t) 
during the transition from the relativistic regime to the non- 
relativistic regime, with initial conditions T s = T) S> m g and 
fj, s = 0 at a = a;. 


Here, N = J dp Anp 2 / (2 tt) 3 x f(p,t) is a normalization 
factor. Besides the conditions of kinetic equilibrium dis¬ 
cussed above, we also need to take into account that ster¬ 
ile neutrino self-interactions freeze out at a time t dcc and 
sterile sector temperature T s = T s dec , after which kinetic 
equilibrium is lost and sterile neutrino momenta are sim¬ 
ply redshifted as a~ 1 (t). This implies for the momentum 
distribution function: 


f(p,t) 


exp [rfftT (v / p 5 +mj-/i >(*))] + 1 

_ 1 _ 

6XP [ TfdA (\J *2 +™s~»s (« dBC )) ] +1 


for t < t dec 
for t > t dec ■ 


(A9) 


Here, p s (t dcc ) is the chemical potential at the time of 
decoupling. We have checked that the exact decoupling 
time only slightly changes the evolution of P, so we re¬ 
gard our solution in Fig. [5] as universal for all parameter 
values of interest. In SecTjlV[ we have for simplicity as¬ 


sumed a sudden decoupling of self-interactions though. 
For the value e 2 s /M 2 ~ 1 MeV -2 chosen there, this leads 
to T, ~ 0.0024 eV, corresponding to a photon tem¬ 
perature of 0.038 eV. 
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